Ultracold RbSr molecules can be formed by magnetoassociation 
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We investigate the interactions between ultracold alkali metal atoms and closed-shell atoms using 
electronic structure calculations on the prototype system Rb+Sr. There are molecular bound states 
that can be tuned across atomic thresholds with magnetic field, and there are previously neglected 
terms in the collision Hamiltonian that can produce zero-energy Feshbach resonances with significant 
widths. The largest eflFect comes from the interaction-induced variation of the Rb hyperfine coupling. 
The resonances may be used to form paramagnetic polar molecules if the magnetic field can be 
controlled precisely enough. 

PACS numbers: 34.50.Cx, 37.10.Pq, 67.85.-d 



There have been enormous advances in ultracold 
molecule formation in the last few years. Alkali metal 
dimers have been formed both by photoassociation and 
by magnetic tuning across zero-energy Feshbach reso- 
nances (magnetoassociation). In both approaches the 
molecules are initially formed in very high-lying vibra- 
tional states, but methods are novi^ emerging to transfer 
the molecules to low- lying vibrational states For 
KRb and Cs2, molecules in the rovibronic ground state 
have been formed coherently from ultracold atoms by 
magnetoassociation followed by stimulated Raman adia- 
batic passage (STIRAP) The availability of sam- 

ples of polar ultracold molecules opens up new possibil- 
ities for exploring polar quantum gases, for developing 
quantum simulators, and for studies of controlled ultra- 
cold chemistry pj| . 

Alkali metal dimers have singlet ground states, which 
cannot be tuned magnetically except through the very 
small magnetic moments of the constituent nuclei. Al- 
though the molecules can be formed in excited triplet 
states, the triplet dimers are likely to be subject to trap 
loss due to fast inelastic and reactive collisions. There 
is therefore great interest in forming ultracold molecules 
that have ground states with unpaired electron spin. 
Molecules that have both electric and magnetic dipoles 
are particularly interesting, because they may be used to 
create topologically ordered states and may have novel 
applications in quantum information storage pT|. 

A very promising class of molecules are those formed 
from an alkali metal atom, with a ground state, and 
a closed-shell species such as an alkaline earth atom or a 
Yb atom, with a ground state. However, the collision 
Hamiltonian that is usually used for ultracold atom colli- 
sions docs not couple the atomic and molecular states in 
such systems, and it is commonly believed that they will 
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not exhibit magnetically tunable Feshbach resonances 
that can be used for ultracold molecule formation. The 
purpose of the present paper is to show that this is not in 
fact true: there are terms in the Hamiltonian that have 
previously been neglected, which can give rise to mag- 
netically tunable Feshbach resonances. If precise enough 
control of magnetic fields can be achieved, these reso- 
nances may be used for molecule formation in these im- 
portant systems. 

We consider the prototype system RbSr. This is par- 
ticularly topical because Bose-Einstein condensation has 
recently been achieved for ^^Sr and ^^Sr ^ and 

Fermi degeneracy for ^''Sr [iH]. Both ®'^Rb and ^^Rb are 
readily condensed. However, our general conclusions ap- 
ply to any system made up of an alkali metal atom and 
a closed-shell species. 

The collision Hamiltonian for a pair of atoms a and b 
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+ Ha + Hb + V{r), (1) 



where r is the internuclear distance, is the angular 
momentum operator for mechanical rotation of the atoms 
about one another. Ha and Hb are the Hamiltonians for 
the free atoms (in an applied field) and V{r) is the in- 
teraction operator. For collision of a pair of alkali metal 
atoms. Ha and Hi, are 
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where is the hyperfine coupling constant for atom a, 
Sa and ia are the corresponding electron and nuclear spin 
operators, and g'™'^ are the ^-factors, and B is the 
magnetic field, whose direction defines the z-axis. The 
interaction operator is usually represented 
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where for a pair of alkah metal atoms Vb(r) = (0|V'(r)|0) 
and Vi{r) = (l|V'(r)|l) are the potential curves for the 
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singlet and triplet states and V'^{r) is a spin-spin term 
that represents the magnetic dipole interaction between 
the electron spins on the two atoms (and may incorporate 
short-range terms due to second-order spin-orbit interac- 
tion) . 

The molecular wavefunctions are conveniently ex- 
panded in an uncoupled basis set \samsa)\ia'mia)\sbmsb) 
\ib'm'ib)\LML). The Hamiltonian is diagonal in the to- 
tal projection quantum number Mtot = Mp + Ml , where 
Mp = msa+fTiia+msb-'rfnib. The only term in the Hamil- 
tonian that has matrix elements ofF-diagonal in L is the 
spin-spin term. However, for a pair of alkali-metal atoms, 
the singlet potential is substantially different from the 
triplet potential (typically a factor of 10 to 20 deeper), 
and the difference produces strong couplings between 
states with the same L and the same Mp. The resulting 
molecular states typically have magnetic moments dif- 
ferent from the atomic states, and may cross them as a 
function of magnetic field. The molecular states are cou- 
pled to the atomic states by either Vi{r) — Vo{r) (when 
the atomic and molecular states have the same values 
of L and Mp) or V'^{r) (when the L or Mp values are 
different). This produces magnetically tunable Feshbach 
resonances that may be used for molecule formation. 

The situation is different when one of the atoms is in a 
state. When Sb = 0, there is no spin-spin interaction. 
There is also only one interaction potential, correspond- 
ing to S* = Sq (a doublet, with S = \, for RbSr). The 
quantum number F = = ± i is well-defined only 
at zero field, but even at finite field the molecular states 
have almost exactly the same mixture of F values as the 
atomic states. For the stable isotopes ^''Sr, ®^Sr and ®®Sr, 
which all have zero nuclear spin, the molecular states 
have almost exactly the same magnetic moment as their 
parent atomic states, so are closely parallel to them as a 
function of magnetic field. Even for ®^Sr, with ib = 9/2, 
the nuclear Zeeman effect is a small perturbation. 

If the interaction operator was really represented by 
Eq. ([3]), with only S — \ and no F^(r), there would 
be no coupling between the atomic and molecular states. 
Although there would still be crossings between atomic 
states and molecular bound states as a function of mag- 
netic field, there would be no coupling between them; the 
resulting Feshbach resonances would have zero width and 
it would be impossible to tune across them adiabatically, 
as required for molecule formation. However, Eq. ^ is 
in reality an approximation, and there are several addi- 
tional ways in which the colliding species interact with 
one another. For RbSr, these additional terms have a 
profound effect. 

The most important additional interaction term comes 
from the fact that the Rb hyperfine coupling constant C, 
is modified when another atom is nearby. We may write 
C(r) = ^Rb + AC(?'). The term A(^(r)«a -Sa is most appro- 
priately considered to be part of the interaction operator 
V{r). In addition, there are short-range contributions 
to V{r) from smaller terms such as (i) the interaction 
eQq{r) between the nuclear quadrupole moment of Rb 
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FIG. 1: Interaction potential V{r) (top panel) with dipole 
moment /i(r) (middle panel) and hyperfine coupling constant 
("(r) (bottom panel). The binding energy and equilibrium dis- 
tance are calculated to be 1000 cm~^ and 4.67 A, respectively. 

and the field gradient at the Rb nucleus caused by the 
distortion of the electron density by Sr; (ii) the dipolar in- 
teractions between the electron spin and the Rb nuclear 
spin; (iii) the spin-rotation interactions between L and 
the Rb electron and nuclear spins. However, the mod- 
ification of is by far the largest effect, as discussed 
below. 

We have carried out high-level electronic structure cal- 
culations of the interaction potential and dipole moment 
of RbSr and the modification of by Sr. The re- 
sults are summarized in Fig. [T] The short-range po- 
tential (from 3 to 10 A) was calculated using the spin- 
restricted coupled-cluster method with single, double and 
approximate triple excitations [CCSD(T)], with the rel- 
ativistic small-core ECP28MDF effective co re p otentials 
and uncontracted basis sets of Lim et al. |l6l . aug- 
mented by 3s3p2(i functions at the bond midpoint. For 
i? > 10 A, we extrapolate using the long-range form 
V(r) ~ —Cer"^ — Csr~^, with semiempirical coefficients 
[181] Cg = 3.762 X 103 Ei^af, and Cg = 4.62 x lO^E^al 
The short- and long-range parts of the potential were 
smoothly connected with the switching function used by 
Janssen et al. [l^. The quantities C(r) and eQq{r) were 
calculated with the relativistic density-functional theory 
(DFT) approach [20] implemented in the ADF program 
using the PBEO functional The asymptotic 

value of ^Rb was underestimated by 6% in the DFT cal- 
culations, so we scaled ^(r) to reproduce the experimen- 
tal atomic value. The quantity A({r) was fitted to the 
Gaussian form (^qc"'^^^'"^''^ , giving parameters a = 0.23 
A-2 and ra = 4.06 A, with Co = -687 MHz for ^''Rh 
and -229 MHz for ^^Rb. This corresponds to a 20% 
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maximum reduction in We also evaluated the dipole 
moment ^{r) using a variety of approaches. The most 
reliable is a finite- field CCSD(T) calculation with the ap- 
proach described above, which gives /i = 1.36 D at the 
equilibrium distance — 4.67 A. However, the precise 
value is sensitive to the level of correlation treatment. 

We have also investigated the smaller additional cou- 
plings described above, in order to verify that they are 
much less important than A({r). All these can cou- 
ple states with AMp ^ when L > 0. The nuclear 
quadrupole coupling and the dipolar coupling between 
Sa and ia can also couple channels with AL ^ and 
can thus mediate Feshbach resonances in s-wave scatter- 
ing for bound states with L 0. Our DFT calculations 
give values for the nuclear quadrupole coupling coupling 
constant eQq{re) = 8 MHz and 3.7 MHz for ^^RbSr and 
^^RbSr respectively, reducing rapidly to zero as the RbSr 
distance increases beyond Tg. The coupling constant for 
the dipolar interaction between Sa and ia is of the order 
of 1 MHz near Vf,. The electronic spin- rotation coupling 
constant can be estimated in terms of the rotational con- 
stant h and the anisotropy (?aniso of the electronic factor 
as 26(7aniso @| , which is approximately 20 MHz near r^- 
The nuclear spin-rotation interaction will be about a fac- 
tor of 10'^ smaller because of the ratio of the nuclear and 
Bohr magnetons. A 20 MHz coupling is potentially sig- 
nificant, but neither the electronic nor the nuclear spin- 
rotation interaction has matrix elements that affect s- 
wave scattering. All these smaller coupling terms, and 
couplings involving the nuclear spin of ®^Sr, are neglected 
in the scattering calculations described below. 

Since the Mi?-changing terms in the collision Hamil- 
tonian are so small, crossings between thresholds and 
bound states with different Mp will not produce Fesh- 
bach resonances wide enough to be measured with cur- 
rent experimental methods. However, if the Rb atom 
is initially in a state that correlates at zero field with 
F = ia — \, then the threshold is crossed by bound states 
with the same Mp but correlating with F = ia + \- The 
operator AC,{r)ia ■ Sa is not diagonal in the field-dressed 
atomic eigenstates, so it can produce Feshbach reso- 
nances at these crossings. The pattern of bound states 
and the crossings that produce Feshbach resonances are 
shown in Fig.[2]for *^Rb and ^''Rb with all the stable iso- 
topes of Sr. Since the hyperfine splittings are a few GHz, 
the bound states responsible for the crossings are bound 
by energies of a few GHz. For RbSr these are states with 
vibrational quantum numbers v — —3, —4 (relative to 
threshold) for magnetic fields up to 0.5 T. 

We next investigated the widths of the resonances pro- 
duced in this way. To locate them, we first carried out 
bound-state calculations as a function of magnetic field 
using the BOUND package 2J] to determine precisely 
the magnetic field at which the crossing occurs. We then 
used the MOLSCAT program [H], modified to handle 
collisions of atoms in magnetic fields [l^, to calculate 
the scattering length a{B) as a function of magnetic 
field near the crossing at a near-zero collision energy 
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FIG. 2; (Color online) Molecular levels (colored) that cross 
atomic thresholds (black) for \Mf\ < ia — ^ states of *^RbSr 
(upper panel) and *^RbSr (lower panel) as a function of mag- 
netic field for different Sr isotopes. The dotted, dashed and 
solid colored lines correspond to v = —2,— 3,— 4 vibrational 
states, respectively. Positions of resonances are marked with 
circles. 



(10~* K). This was then fitted to the functional form 
a{B) ~ flbg [1 + A/(i? — -Brcs)] to extract the resonance 
position Bros, the width A, and the background scatter- 
ing length flbg. 

The parameters of selected Feshbach resonances are 
given in Table HI It should be noted that our calcu- 
lated interaction potential is not accurate enough to pre- 
dict the positions of the highest bound states (or the 
scattering length) correctly for a specific isotope. How- 
ever, there are enough different isotopes of Sr available 
that one of them is likely to display crossings of each of 
the types shown in Fig. [2] A measurement of the scat- 
tering length or the binding energy of one of the near- 
dissociation states for any isotopic species will allow reli- 
able calculations of the resonance positions for all isotopic 
combinations. 

States with Mp < (except Mp = —ia — have 
energy maxima or minima at -Bturn = ^QMp / [g^^i^)^ 

1/2 

where they are separated hy d = C, \{ia + ^)^ — 

For levels with binding energies in the range d < 
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TABLE I: Calculated properties of RbSr Feshbach resonances 
at F — ia — ^ thresholds arising from bound states supported 
by the F — ia + ^ thresholds. 
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Ev < (*a + ^)C: there are two crossings with the lower 
threshold: one between and Bturn and the second be- 



tween -Bturn and 2i?turn- Convcrsely, for levels with bind- 
ing energies \Ey \ > (zq -|- i)^, there is only one crossing, 
at a field B > 2i?turn- For levels with Mp > 0, there is 
only one crossing for each vibrational state with binding 
energy larger than (zq -|- l/2)(. 

The resonance width is proportional to the square 
of a bound-continuum matrix element {v\A^{r)ia ■ 
Sa I continuum). The radial part of this is proportional to 
the amplitude of the bound-state wavefunction at short 
range, which varies as where n = 6 is the 

power of the leading term in the long-range potential. 
Since A({r) is itself roughly proportional to C, the widths 
of the low-field resonances may generally be expected 
to increase with ( approximately as C®/"^. Because of 
this, the widths observed for *^Rb in Table U are gen- 
erally larger than those for ®^Rb. However, there are 
other factors involved such as the likelihood of obtain- 
ing low-field resonances at all (which decreases with ^) 
and the magnitude of the off-diagonal matrix element 
of ia ■ Sa, which increases linearly with B at low field 
and then levels off above iJturn to a value proportional 

1/2 

to [{ia + 5)^ — Mp] . Because of this, the lowest-field 
resonances in Table U are all very narrow. Resonance 
widths are also enhanced in cases where the background 
scattering length is large. 

The widest and most promising resonances for 
molecule production are of two types. The first are those 
that occur for Mp < states that cross the lower thresh- 
old twice, as occurs for v — —3 for *^Rb®''Sr and v = —4 
for ®''Rb®'*Sr in Fig. [2] The higher-field resonance is al- 
ways the wider of the pair because of the ia ■ Sa matrix 
element discussed above. The second are those for which 
the background scattering length is large, as for ^^Rb^^Sr 
and ^^Rb^'^Sr on the current potential. These two effects 
combine for the 1959 G resonance for *^Rb®^Sr in Table 
Uto give a width as high as 122 mG. 

In conclusion, we have investigated the interactions be- 
tween Rb and Sr atoms and have identified a new mech- 
anism that can produce magnetically tunable Feshbach 
resonances in collisions of ultracold molecules. These 
Feshbach resonances could be used to produce ultracold 
molecules that would have both electric and magnetic 
dipole moments in their ground states. The resonances 
arise from the modification of the Rb hyperfine coupling 
by the presence of another atom. The mechanism is gen- 
eral and may produce magnetically tunable Feshbach res- 
onances in any system in which an atom with electron 
spin collides with a closed-shell atom. 

This work is supported by EPSRC under collaborative 
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CORES Programme EuroQUAM. 



[1] J. M. Sage, S. Sainis, T. Bergeman, and D. DeMille, [2] K.-K. Ni, S. Ospelkaus, M. H. G. de Miranda, A. Pe'er, 
Phys. Rev. Lett. 94, 203001 (2005). B. Neyenhuis, J. J. Zirbel, S. Kotochigova, P. S. Julienne, 



5 



D. S. Jin, and J. Yc, Science 322, 231 (2008). 
[3] F. Lang, K. Winkler, C. Strauss, R. Grimm, and 
J. Hockcr Denschlag, Phys. Rev. Lett. 101, 133005 
(2008). 

[4] J. G. Danzl, E. Haller, M. Gustavsson, M. J. Mark, 

R. Hart, N. Bouloufa, O. Dulieu, H. Ritsch, and H.-C. 

Nagerl, Science 321, 1062 (2008). 
[5] M. Vitcau, A. Chotia, M. AUogrini, N. Bouloufa, 

O. Dulieu, D. Comparat, and P. Pillet, Science 321, 232 

(2008). 

[6] J. Deiglmayr, A. Grochola, M. Repp, K. Mortlbauer, 
C. Gliick, J. Lange, O. Dulieu, R. Wester, and M. Wei- 
demiiller, Phys. Rev. Lett. 101, 133004 (2008). 

[7] M. J. Mark, J. G. Danzl, E. Haller, M. Gustavs- 
son, N. Bouloufa, O. Dulieu, H. Salami, T. Bcrgeman, 
H. Ritsch, R. Hart, et al., Appl. Phys. B 95, 219 (2009). 

[8] C. Haimberger, J. Kleinert, P. Zabawa, A. Wakin, and 
N. P. Bigelow, New J. Phys. 11, 055042 (2009). 

[9] J. G. Danzl, M. J. Mark, E. Haller, M. Gustavsson, 
R. Hart, J. Aldegundo, J. M. Hutson, and H.-C. Nagerl, 
Nature Physics 6, 265 (2010). 
[10] L. D. Carr, D. DcMillo, R. V. Krems, and J. Ye, New J. 

Phys. 11, 055049 (2009). 
[11] A. Micheli, G. K. Brennen, and P. ZoUer, Nature Physics 
2, 341 (2006). 

[12] S. Stellmer, M. K. Tey, B. Huang, R. Grimm, and 
F. Schrcck, Phys. Rev. Lett. 103, 200401 (2009). 

[13] Y. N. Martinez de Escobar, P. G. Mickelson, M. Yan, 
B. J. DcSalvo, S. B. Nagel, and T. C. Killian, Phys. Rev. 
Lett. 103, 200402 (2009). 

[14] P. G. Mickelson, Y. N. Martinez de Escobar, M. Yan, 



B. J. DoSalvo, and T. C. Killian, arXiv:1003.3867 (2010). 
[15] B. J. DeSalvo, M. Yan, P. G. Mickelson, Y. N. Mar- 
tinez de Escobar, and T. C. Killian, arXiv: 1005.0668 
(2010). 

[16] I. S. Lim, H. StoU, and P. Schwerdtfeger, J. Chem. Phys. 

124, 034107 (2006). 
[17] 1. S. Lim, P. Schwerdtfeger, B. Metz, and H. StoU, J. 

Chem. Phys. 122, 104103 (2005). 
[18] J. Mitroy and M. W. J. Bromley, Phys. Rev. A 68, 

052714 (2003). 

[19] L. M. C. Janssen, G. C. Groenenboom, A. van der 
Avoird, P. S. Zuchowski, and R. Podeszwa, J. Chem. 
Phys. 131, 224314 (2009). 

[20] E. van Lenthe, E. J. Baerends, and J. G. Snijders, J. 
Chem. Phys. 101, 9783 (1994). 

[21] G. to Velde, F. M. Bickclhaupt, S. J. A. van Gisbergen, 

C. Fonseca Gucrra, E. J. Baerends, J. G. Snijders, and 
T. Ziegler, J. Comput. Chem. 22, 931 (2001). 

[22] C. Adamo and V. Barone, J. Chem. Phys. 110, 6158 
(1999). 

[23] R. F. Curl Jr., Mol. Phys. 9, 585 (1965). 

[24] J. M. Hutson, E. Ticsinga, and P. S. Julienne, Phys. Rev. 
A 78, 052703 (2008). 

[25] J. M. Hutson and S. Green, MOLSCAT computer pro- 
gram, version I4, distributed by Collaborative Computa- 
tional Project No. 6 of the UK Engineering and Physical 
Sciences Research Council (1994). 

[26] M. L. Gonzalez-Martmez and J. M. Hutson, Phys. Rev. 
A 75, 022702 (2007). 



